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Introduction 


This  report  constitutes  the  final  report  on  work  performed 
under  Air  Force  Office  of  Scientific  Research  Contract  No.  F49620-81-C-0004 
(Amendment  P00001),  Project  No.  2301/A7  entitled  "Theoretical  Investiga¬ 
tion  of  Thick-Target  Nuclear  Signatures". 

The  subject  of  the  study  concerns  predictions  of  y-rays  and 
neutrons  emitted  In  the  backward  directions  from  thick-targets  of  various 
material  compositions  when  bombarded  by  various  types  of  ions  (hydrogen 
through  lithium)  In  the  energy  range  from  about  50  MeV  to  several  hundred 
MeV. 

The  original  scope  of  the  study  consisted  of  the  following 

tasks : 

1.  extend  the  capabilities  of  previous  nuclear  interaction 
models  and  radiation  transport  computer  codes  to  treat 
essentially  arbitrary  target  composition  and  to  allow 
heavier  Ion  beam  masses  (through  lithium  Ions), 

2.  carry  out  calculations  of  the  "backscattered"  neutron 
and  y-ray  spectra  for  specific  cases  for  which  experi- 
ments  (under  a  separate  program)  were  to  be  performed, 

3.  perform  systematic  theoretical /experimental  comparisons 

of  yields  and  spectra  to  access  the  accuracy  of  the  models 
and  calculatlonal  methods,  and 

4.  assess  the  detectability  of  the  y-  and  neutron- radiation, 
including  background  radiation  effects. 

(The  detailed  statement-of-work  for  the  original  study  is  included  here 
In  the  appendix.) 

Due  to  unforeseen  circumstances,  the  planned  associated 
experimental  program  was  not  carried  out,  so  the  performance  of  Tasks  2-4 
above  were  not  possible.  Thus,  only  part  of  Task  1  above  was  carried  out 
under  the  present  study.  This  work  consisted  of  Initial  nuclear  model 
modifications  and  extensions  to  Include  fast-fission  reactions  for  uranium 
targets,  as  described  In  Section  2. 
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Subsequent  to  this  study,  under  separate  support  from  the  Air 
Force  Weapons  Laboratory  (AFVIL),  work  under  Task  1  was  completed,  the 
computer  codes  were  made  operable  on  computers  at  the  AFWL,  and  calcu¬ 
lations  were  carried  out  for  several  beam  parameters.  A  few  of  these 
results  are  Included  here  as  Section  3  to  Indicate  that  the  predicted 
y  and  neutron  yields  and  spectra  appear  to  be  reasonable  qualitatively, 
but  to  date  there  has  not  been  experimental  verification. 


2.  Nuclear  Model  Extensions 

The  basic  method  employed  for  calculating  the  neutron  and  y-ray 
yields  from  high-energy  Ion  bombardnent  of  thick  targets  Is  to  use  theoreti¬ 
cal  nuclear  Interaction  models  and  thick  target  radiation  transport  com¬ 
puter  codes.  Monte  Carlo  methods  are  used  In  computing  both  the  nuclear 
Interaction  products  and  the  radiation  transport. 

The  work  here  Is  an  extension  of  previous  methods  that  have  been 
developed  for  lower-mass  Ions,  as  summarized  In  Section  2.1  below.  Modi¬ 
fications  made  In  this  study  to  allow  Incident  Ions  up  to  lithium  are 
described  i*  Section  2.2,  and  modifications  to  incorporate  fast  fission 
reactions  to  allow  predictions  for  uranium  targets  are  described  In  Section 
2.3. 

2.1  Background 

Consider  an  energetic  Ion  beam  having  a  kinetic  energy  of,  say, 
a  few  hundred  MeV  and  low-mass  (In  the  range  from  *H  to  7L1)  bombarding 
a  target  where  the  target  thickness  Is  £  the  Ion  range.  The  main  physical 
processes  to  be  simulated  are  as  follows: 

1.  Inside  the  target,  the  primary  beam  loses  energy  by 
atomic  processes  (ionization  and  excitation)  and  under¬ 
goes  nuclear  Interactions  ("spallation"  collisions). 

2.  These  spallation  reactions  prodtxe:  (a)  high-energy 
neutrons  and  protons  (up  to  the  energy /nucleon  of  the  ion), 
which  travel  to  produce  additional  spallation  reactions, 
inducing  a  "hadron  cascade"  Inside  the  target  material, 

(b)  low-energy  particles  (deuterons,  tritons,  etc.)  which 
have  very  short  ranges,  and  are  not  of  Interest  here,  and 

(c)  "spallation"  y-rays  from  nuclear  de-excitation. 

3.  The  low-energy  neutrons  (^  15  MeV)  produced  can  also  undergo 
Inelastic  scattering  (n,  n*y)  and  capture  (n,y)  reactions 

to  produce  an  Additional  source  of  y-rays. 

The  main  concern  here  is  to  simulate  In  detail  these  Interaction  and  trans¬ 
port  processes  to  determine  the  energy  spectrum  of  those  neutrons  and  y-rays 
which  escape  from  the  beam-entrance  face  of  the  target. 


Nuclear  Models 

We  divide  nuclear  collisions  into  two  categories:  (a)  spallation 
collisions,  which  we  take  to  be  collisions  induced  by  all  particles  (neu¬ 
trons,  protons,  and  primary  ions)  having  energies  £  20  MeV  per  nucleon, 
and  (b)  low-energy  (£  20  MeV)  neutron  collisions.  The  reason  for  this 
division  Is  that  experimental  data  are  available  for  treating  low-energy 
neutron  reactions  (from  extensive  research  related  to  nuclear  reactors),  so 
we  can  utilize  the  rather  extensive  experimentally  derived  cross  section 
libraries  In  this  energy  range.  For  the  spallation  reactions,  adequate 
experimental  data  are  not  available,  so  we  resort  a  theoretical  model,  as 
described  below.  For  charged  particles  below  ^  20  MeV/nucleon,  It  Is  assumed 
that  the  particles  stop  and  come  to  rest  before  undergoing  nuclear  collision. 

The  Intranuclear-cascade-evaporation  (ICE)  model  has  been  applied 
successfully  for  some  years  for  treating  spallation  reactions  Induced  by' . 
nucleons  and  charged  plons  (e.g..  Ref.  1  ).  This  model  assumes  that  the 
reaction  takes  place  In  two  phases.  In  the  first  phase  (Intranuclear  cascade 
step)  the  collisions  which  take  place  Inside  the  nucleus  between  the  Inci¬ 
dent  particle  the  Individual  neutrons  and  protons  Inside  the  nucleus  are 
computed  using  Monte  Carlo  methods,  taking  Into  account  the  Fermi  motion 
of  the  nucleons  and  forbidden  collisions  due  to  the  Pauli  principle.  An 
Important  advantage  of  this  model  Is  that  the  principal  Input  data  re¬ 
quired  are  particle-particle  cross  sections  (nucleon-nucleon  and  pion- 
nucleon)  which  are  rather  well  known  experimentally.  The  residual  ex¬ 
citation  energy  after  the  Intranuclear  cascade  Is  assumed  to  be  uniformly 
distributed,  and  a  statistical  ("evaporation")  model  Is  then  applied  as 
a  second  step  to  determine  further  particle  emission  via  de-excitation. 

This  model  has  also  been  applied  to  treat  heavier  Incident  particles 
(deuterons  through  alpha  particles).  Ref.  2,  and  was  extended  In  this  study 
to  treat  projectile  Ions  through  lithium  (Section  2.2). 

The  "standard"  ICE  model  assumes  that  the  residual  excitation 
energy  left  after  particle  emission  Is  no  longer  energetically  possible  by 
evaporation  Is  dissipated  by  y-ray  emission,  but  the  y-ray  spectrum  Is  not 
determined.  Extensions  of  the  model  have  been  made  previously  to  determine 
the  energy  spectrum  of  these  spallation  y-rays  by  computing  the  transition 
probabilities  and  energies  between  energy  levels  during  de-excitation 
(Ref.  3  ). 


Radiation  Transport 

For  the  transport  of  the  Incident  Ions  and  high-energy  nucleons, 
a  modified  version  of  the  Monte  Carlo  radiation  transport  code  HETC  (Ref. 

4  )  Is  used.  This  code  Incorporates  the  ICE  spallation  model,  and  com¬ 

putes  the  trajectories  and  collisions  for  all  secondary  particle  production. 

The  high-energy  transport  code  also  computes  the  low-energy  neu¬ 
trons  produced  In  spallation  collisions.  The  energy,  direction,  and 
spatial  point  of  each  low-energy  neutron  (ti  20  MeV)  produced  In  spallation 
collisions  during  the  high-energy  radiation  transport  Is  written  on  an 
output  tape.  This  tape  Is  then  used  as  an  Input  neutron  source  to  the  low- 
energy  transport  code  MORSE  (Ref.  5  } .  The  MORSE  code  Is  then  used  to 
perform  a  coupled  neutron/y-ray  transport  calculation  -  l.e.,  the  y-rays 
produced  by  low-energy  neutron  Interactions  are  also  transported.  The 
nuclear  Interaction  cross  sections  needed  for  the  MORSE  calculations  can. . 
be  taken  from  several  available  compilations  which  are  maintained  and  rou¬ 
tinely  updated  with  newly  measured  data,  such  as  the  Evaluated  Nuclear  Data 
File  (Ref.  6)... 

The  high-energy  transport  calculation  also  provides  the  residual 
excitation  energy  and  type  of  nucleus  (charge  and  mass  numbers)  left  after 
spallation  reactions.  This  Information  for  each  collision  Is  used  to  com¬ 
pute,  in  a  separate  computer  program,  the  spectrum  of  y-rays  produced. 

These  spallation  y-rays  are  then  used  as  a  source  for  the  MORSE  code  to 
determine  the  transport  of  these  y-rays  throughout  the  target. 

Both  the  high-energy  and  low-energy  transport  codes  are  general 
purpose,  providing  a  "record"  of  all  essentially  events  which  take  place 
during  the  Monte  Carlo  simulation.  Analysis  programs  of  these  recorded 
histories  are  then  written  to  extract  the  relevant  Information  for  a  par¬ 
ticular  problem  of  Interest.  For  present  interests,  the  analysis  programs 
are  written  to  compute  the  spectra  (and  Integral)  of  neutrons  and  y-rays 
escaping  the  target  through  the  beam-entrance  surface. 

A  summary  of  the  calculatlonal  steps  Involved  Is  given  In  Figure  1. 
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I  Input:  •  Bean  Parameters 


e  Target  Information  | 


2.2  Intranuclear-Casecade-Evaporation  Model  Modifications 

For  Ion-Nucleus  Collisions 

In  extending  the  ICE  model  to  predict  the  products  from  reac¬ 
tions  Induced  by  projectiles  of  A  >  1,  the  basic  assunption  is  made  that 
the  projectile  can  be  considered  as  a  cluster  of  independent  nucleons. 

The  reaction  is  then  computed,  In  effect,  as  the  sum  of  the  effects  of 
these  Independent  nucleons  striking  the  nucleus.  This  approximation  would 
appear  to  be  quite  valid  for  Ions  only  slightly  higher  than  nucleons. 

However,  this  approximation  Is  expected  to  become  progressively  worse  as 
the  projectile  Ion  mass  Increases,  where  cascade  reactions  can  be  induced 
in  the  projectile  particle  as  well  as  the  target  nucleus.  In  the  present 
work  we  have  made  modifications  to  the  ICE  model  using  the  Independent 
nucleon  approach  to  allow  projectiles  up  to  A  «  10,  although  the  heaviest 
ion  for  which  complete  calculations  have  been  made  to  date  using  the  method 
outlined  in  Figure  1  has  been  for  7Li. 

To  determine  the  collision  site  and  the  type  of  struck  nucleus  in 
the  target  material  composition,  an  empirical  formula  for  the  nonelastic 
nuclear  cross  section  is  used:  a  *  ir(R^  +  Rj  -  6)2,  where  Rj  *  roAj^3  is 
the  radius  of  the  target  nucleus,  R^  *  2.19  A-f^3  x  10"13  cm  is  the  ion  radius, 
r0  -  1.42  x  10~13  cm,  and  6  =  0.80  x  10*13  cm. 

Next,  the  impact  parameter  of  the  center-of-mass  of  the  ion  is 
determined  by  selecting  a  point,  with  uniform  probability,  over  a  projected 
circular  area  of  radius  R-j  +  Rj.  To  determine  the  spatial  coordinates  at 
which  each  nucleon  of  the  projectile  strikes  the  target  nucleus,  a  con¬ 
figuration  for  the  cluster  relative  to  the  projectile  conter-of-mass  must 
be  assumed.  We  assume  that  one-half  of  the  nucleons  are  located  at  each 
end  of  a  line  a  distance  R^  from  the  projectile  center  of  mass,  with  proba¬ 
bility  of  0.5  that  the  additional  nucleon  for  odd  Is  located  at  either  end. 
The  projectile  cluster  Is  allowed  to  have  a  uniform  azimuthal  orientation. 
Thus,  this  procedure  defines  the  position  on  the  surface  of  the  target 
nucleus  which  each  nucleon  of  the  projectile  enters.  The  kinetic  energy 
of  each  nucleon  as  It  enters  the  nucleus  is  taken  to  be  (E  -  Eg)/A|,  where 
E  and  Eg  are  the  kinetic  and  binding  energies  of  the  ion. 

An  Intranuclear  cascade  calculation  is  then  performed  for  each 
projectile  nucleon,  and  the  sum  of  the  produced  secondary  particles  represents 
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the  cascade  products  for  the  Ion  Interaction.  The  A,  Z,  and  E*  (excitation 
energy)  after  the  cascade  phase  Is  then  used  as  input  for  an  evaporation 
calculation  in  the  same  manner  as  In  the  standard  model  for  nucleon- 
nucleus  collisions. 

For  targets  containing  hydrogen,  ion  collisions  with  proton 
nuclei  are  treated  using  the  same  nucleon-nucleon  cross  sections  as  are 
used  in  the  Intranuclear  cascade  portion  of  the  calculation  together  with 
the  algorithm  given  In  Ref.  Z.  The  projectile  ion  energy  loss  by  atomic 
ionization  and  excitation  is  determined  by  scaling  the  stopping  values  for 
protons,  as  described  in  Ref.  2. 

2.3  Fast  Fission  Mechanism 

For  particle  Interactions  with  very  heavy  target  nuclei  (e.g. , 
uranium),  the  two-step  intranuclear  cascade  plus  evaporation  model  dis- 
cussed  above  should,  according  to  experimental  evidence,  be  modified  to 
take  Into  account  an  additional  de-excitation  mechanism:  "fast"  fission. 
That  is,  after  the  Intranuclear  cascade  step,  there  Is  "competition"  at 
each  stage  of  de-excitation  (l.e.,  after  each  particle  emission)  as  to 
whether  the  residual  nucleus  continues  to  evaporate  particles  from  a  single 
fragment  or  whether  the  residual  nucleus  fissions  Into  two  fragments,  with 
each  fission  fragment  continuing  de-excitation  by  evaporating  particles. 

The  overall  emitted  particle  multiplicities  and  spectra  from  the  Inter¬ 
action  are  expected  to  be  somewhat  different  If  fission  takes  place. 

The  probability  that  fission  will  occur  given  a  nonelastic  Inter¬ 
action  depends  upon  the  particle  type,  particle  energy,  and  target  nucleus. 
For  example,  for  protons  of  ^  100  MeV,  the  fission  probability  is  %  0.8 
for  U2*8,  Jfc  0.5  for  Th2*2,  and  ^  0.05  for  Bi208,  with  the  fission  proba¬ 
bility  varying  roughly  as  Z2/A  of  the  target  nucleus  (Ref.  7).  Thus, 
for  present  Interests  high-energy  fission  is  most  Important  only  for  tar¬ 
get  compositions  containing  uranium.  In  general,  the  probability  of 
fission  Increases  with  particle  mass  and  energy. 

In  order  to  provide  the  computational  capability  of  predicting 
escaping  neutrons  and  ganma-rays  from  thick  targets  containing  uranium, 
the  ICE  model  was  modified  to  incorporate  the  fast  fission  mechanism. 

The  basic  procedure  and  assumptions  follow  those  of  References  8  -  10, 
and  the  major  steps  are  sunmarlzed  below. 
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The  main  Items  to  be  specified  are  the  probability  of  fission 
at  each  step  of  de-excitation  and  a  description  of  the  fission  fragments 
(l.e.,  kinetic  energy,  excitation  energy,  Z,  and  A)  If  fission  takes  place. 
The  fission  probability  Is  written  as 

pf  -  (i  -  yrf)- 

where  rn  and  rf  are  the  widths  for  neutron  emission  by  evaporation  and  for 
fission,  respectively.  In  the  actinide  region  (Z  >  90),  the  empirical  for¬ 
mula  of  Vandenbosh  and  Hulzenga  (Ref.  11)  Is  used, 

r_ 

log  ■  *(Z)A  +  $(Z), 
rf 

where  the  coefficients  $  and  4>  are  based  on  systematlcs  and  have  been  tab¬ 
ulated  In  Ref.  12  .  Pf  Is  assumed  to  be  Independent  of  excitation  energy 
E*  for  E*  >  6  MeV,  and  Pf  *  0  for  E*  <  6  MeV.  For  residual  nuclei  In 

the  subactlvlte  region  (Z  <  90),  statistical  model  fits  to  experimental 
data  are  used: 

Tn  -  0.352[1.68Iq  +  1.93A1/sIi  +  AV»  (0.761!  -  0.05IQ)] 

I,  "  |jat<Sh-1.0).S"  +  1] 

Ii  -  —  [(6.0  -  6.0Sn  +  2.0Sp)  esn  +  Sft  -  6.0] 
an 

Sn  »  2.0  C(an  (E*  -  B)]1* 

B  '■  separation  energy  minus  paring  energy 

an  -  (A  -  1.0)/8 

rf  -  C(Sf  -  1.0)  e5f  +  1.0]/af 

Sf  -  2.0[af(E*  -  Ef) ]** 

Ef  -  B  -  321.2  -  16.7(Z2/A)  +  0.219(Z2/A)2 

af/an  -  1.09  +  0.011 (Z2/A  -  31.1 )2. 

The  basic  equations  above  are  based  on  the  statistical  model  (e.g..  Ref.  13). 
The  parameters  which  have  been  fit  to  experimental  data  are  the  fission  bar¬ 
rier  Ef  and  8f/an  (ratio  of  level  density  parameter  for  fission  to  level 
density  parameter  for  neutron  emission). 


The  kinetic  energy  of  the  fission  fragments  arises  from  mutual 
coulomb  repulsion,  so  the  sum  of  fragment  energies  Ek  Is  expected 
to  correlate  with  the  coulomb  repulsion  parameter  Z2/A  .  From  the  fits 
to  experimental  data  made  In  Ref.  14  , 

Ek(MeV)  -  0.1065(Z2/AVs)  +  20.1 

The  recoil  energy  for  each  fragment  Is  then  given  by  the  relations  Ek  « 
[AZ/(A1  *  Az)]  £k  and  Ekj  ■  [Ai/(A!  *  A2)]  Ek.  The  excitation  energy 
of  the  fragments  Is  then  determined  by  energy  conservation  -  l.e., 

Et  -  M(A,Z)  +  E*  -  M(At.Zi)  -  M(A2,Z2) 

where  M  denotes  nuclear  mass  and  E*  Is  the  excitation  energy  of  the  fission¬ 
ing  nucleus.  The  fragment  excitation  energies  then  are 

E*  -  [ Ai/(A j  +  A2)]  (Et  -  Ek)  and 

Et  -  [A2/(A2  +  A2)3  (Et  -  Ek). 

The  most  appropriate  method  for  specifying  fragment  masses  Is 
not  clear,  either  from  theory  or  measurements.  The  shape  of  fission  mass 
distribution  can  be  either  symmetric  (one  peak)  or  asymmetric  (two  peaks), 
depending  on  the  fissioning  nucleus  and  Its  excitation  energy  (e.g. ,  Ref.  7  ). 
In  the  actinide  region  the  data  of  Ref.  14  Is  utilized,  which  shows  that 
the  mean  mass  of  the  heavy  fragment  Is  essentially  constant  at  about  140 
Independent  of  the  fission  nucleus  mass,  Ap  (for  Ap  %  225  to  %  255),  and 
the  light  fragment  mass  Increases  linearly  from  about  90  to  12C  as  Ap 
Increases.  In  the  subactinide  region,  the  mean  mass  of  each  fragment  Is 
about  Ap/2,  with  a  normal  distribution  having  a  half-wldth-at-half-maxlmum, 
<W>,  of  (Ref.  15):  <W>  -  E*  -  Ef  -  7  (a.m.u.). 

The  procedure  for  treating  high-energy  fission  reactions  then  Is 
as  follows:  First,  the  Intranuclear  cascade  phase  of  the  calculation  Is 
performed  In  the  usual  manner  (l.e.,  this  "fast"  part  of  the  Interaction  Is 
the  same  whether  fission  takes  place  or  not).  The  residual  nucleus  A,Z,  and 
E*  Is  then  used  with  the  probability  of  fission  defined  above  to  "select" 
(In  a  Monte  Carlo  sense)  whether  the  nucleus  at  this  stage  de-excltes  by 
undergoing  fission  or  evaporation.  If  fission  Is  selected,  a  description 
of  the  fragments  Is  determined  using  the  methods  described  above,  and 
evaporation  calculations  are  performed  for  each  fragment  to  compute  par¬ 
ticle  emissions  during  complete  fragment  de-excitation.  If  evaporation  was 
selected,  an  evaporation  calculation  for  the  cascade  residual  nucleus  Is 


performed  for  a  single  particle  emission,  and  the  probability  of  fission  Is 
then  tested  again.  Thus,  the  evaporation-fission  competition  Is  applied 
at  each  step  of  de-excitation  -  l.e.,  after  each  evaporated  particle. 

The  fast- fission  reaction  mechanism  as  described  above  has  been 
programmed  In  the  HETC/LHI  code,  but  Is  largely  wtested  at  this  time. 
(Low-energy  fission  reactions,  due  to  neutrons  below  %  20  MeV  produced  In 
high-energy  collisions,  can  be  taken  Into  account  as  one  of  the  standard 
options  of  the  MORSE  transport  code.) 


3.  Results  From  Transport  Calculations 

While  not  part  of  the  present  study,  the  nuclear  model  develop¬ 
ment  described  In  Section  2  has  been  applied  for  one  case  under  an  AFWL 
sponsored  study  (Ref. 16  ),  and  a  few  of  these  results  are  Included  here  to 
Illustrate  applications  of  the  models. 

The  case  considered  was  L{7  Ions  of  380  and  840  MeV  incident  on 
a  thick  alunlnum  target.  (These  energies  correspond  to  ranges  In  aluminum 
of  about  3  and  10  g/cm2.)  Radiation  transport  calculations  were  carried 
out  using  the  code  system  shown  previously  in  Figure  1.  The  spectra  of 
neutrons  and  gamma-rays  escaping  the  target  In  the  backward  hemisphere 
are  shown  In  Figure  2.  The  "backscattered"  gamma-ray  and  neutron  yields 
for  Lj  Ions  are  compared  in  Figures  3  and  4  with  results  for  other  ions. 

Thus,  the  Li  interaction  model  appears  to  work  properly  based 
on  systematlcs,  but  there  are  no  measured  data  available  for  quantitative 
evaluation. 
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FLUENCE  CM  *  MEV'  (INCIDENT  PARTICLE)" 


Figure  2.  Neutron  and  gamma- ray  spectra  emitted  over  all  backward  directions 
for  340  and  840  MeV  LI  beams  Incident  on  an  aluminum  target 
(from  Ref.  16). 
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Figure  3.  Dependence  of  the  number  o?  "backscr.ttered 
Ion  energy  (from  Ref.  16). 
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Figure  4.  Like  Figure  3  but  for  Neutrons  (from  Ref.  16) 
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a.  The  contractor  will  exercise  his  best  efforts  to  perforin  a  theoretical 
investigation  of  the  nuclear  return  signature  of  thick  targets  bombarded  by 
energetic  light  nuclei.  The  study  will  generally  include  both  calculational  and 
assessment  tasks.  Specific  tasks  to  be  accomplished  will  include  the  following: 

(1)  Implement  modifications  and  improvements  in  the  HETC/LHI  code 
system  required  to  efficiently  perform  parametric  yield  predictions  for  a  range 
of  light  beam  ions,  energies,  and  targets,  and  shall  deliver  an  operating  version 
of  the  modified  HETC/LHI  code  system  source  listing  to  AFWL/NTYP. 

(2)  Perform  a  series  of  specific  calculations  of  neutron  and  gamma  ray 
thick- target  yields.  Representative  light,  intermediate,  and  heavy  target  nuclei 
will  be  considered,  and  beam  particles  of  several  masses  between  W-  and  Li&  will 
be  included.  Energies  of  the  incident  beam  ions  will  be  chosen  to  provide  coverage 
consistent  with  penetration  depths  similar  to  the  proton  energy  region  50-200  MeV. 
Specific  energies  needed  for  validation  by  comparison  with  experiment  will  be 
included. 


(3)  Organize  and  interpret  the  calculational  results.  Comparisons 
with  available  experimental  data  will  be  made,  and  interpretation  of  possible 
discrepancies  and  uncertainties  will  be  made  in  terms  of  their  significance  to 
model  validation. 

(4)  Carry  out  a  preliminary  detectability  assessment,  with  the  theoretical 
results  as  a  basis.  Background  sources  will  be  identified  and  quantified.  Results 
of  the  detectability  assessment  will  include  comparisons  with  estimated  X-ray  yields 
and  backgrounds. 


(5)  Document  the  results  of  the  study  as  appropriate.  Fundamental  results 
will  be  made  available  to  the  technical  ccrmunity  through  the  open  technical  literature 
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